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Abstract 

The behaviour of the space-like string tension in the high temperature phase is 
studied. Data obtained in the Z-2 gauge model in (2+1) dimensions are compared with 
predictions of a simple model of a fluctuating flux tube with finite thickness. It is 
shown that in the high temperature phase contributions coming from the fluctuations 
of the flux tube vanish. As a consequence we also show that in (2+1) dimensional 
gauge theories the thickness of the flux tube coincides with the inverse of the decon- 
finement temperature. 
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1 Introduction 



It is well known that in finite temperature Lattice Gauge Theories (LGT) the "spatial" 
string tension (namely that extracted from Wilson loops orthogonal to the compactified 
imaginary time direction) is no longer the order parameter of confinement and is, in 
general, different from zero even in the deconfined phase. Recently, in a set of interesting 
papers Q [2], || the behaviour of this spatial string tension was studied for SU(2) and 
SU(3) gauge theories in (2+1) and (3+1) dimensions, leading to new ideas and to a better 
understanding of the high temperature behaviour of LGT's. At the same time, it was 
observed in ||] in the case of SU(2) in (2+1) dimensions (and it will be confirmed in the 
present paper, in the case of the gauge Ising model), that the inverse of the deconfinement 
temperature 1/T C almost coincides with the thickness L c of the flux tube joining the 
quark-antiquark pair in the confined phase. This implies that looking to space-like Wilson 
loops in the deconfined phase is almost equivalent to probing the interior of the flux 
tube and gives us a powerful tool to test the effective flux tube models of confinement. 
Moreover, the fact that 1/T C and L c have similar values is a quite interesting phenomenon 
in itself and deserves further investigation. 

In this paper we want to further pursue this analysis, comparing some predictions of 
the effective flux tube picture of LGT's with a set of high statistics Monte Carlo data on 
the 3d Ising gauge model. The main results of our analysis are the following: 

a] In the deconfined region (T > T c ) the space-like string tension increases as the tem- 

perature increases. This trend is rather impressive: almost two orders of magnitude 
are gained moving from the deconfinement point to the highest temperature we can 
measure . Moreover in the deconfined phase two distinct regimes can be identified: 
a first, rather smooth, crossover region (T c < T < 2T C ) and a second, high tem- 
perature regime (T > 2T C ), where the string tension scales with the temperature, 
which is the only remaining dimensional scale in this regime. In particular, this 
second region (T > 2T C ) is the one in which dimensional reduction has been shown 
to apply in the case of the (3+1) dimensional SU(2) gauge theory |J. All these 
results are in complete agreement with the scenarios described in ref.s pi [2L 0], and 
can be understood within the framework of a "compressed" flux tube model. 

b] In the range T c < T < 2T C the simple picture of a compressed flux tube of uniform 

flux density for the space-like Wilson loops is quite accurate. It predicts a string 
tension rising linearly with the temperature and allows a rather precise estimate of 
the thickness L c of the flux tube. 

c] Under the assumption that in the deconfined region (T > T c ) the space-like flux tube 

is "frozen" and the contribution of its transverse quantum fluctuations is zero, it is 
possible to show that L c and 1/T C must coincide. Suitable ratios of Wilson loops 
can be introduced to test this assumption, which turns out to be in good agreement 
with our Monte Carlo simulations. These ratios can be used as order parameters 
for this phase, being zero in the high temperature region and different from zero in 
the confining phase. 
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d] There is an impressive agreement between some dimensionless ratios of physical ob- 
servables that we obtain in the Z2 case and those obtained by Teper ||, |j for 
SU(2). This fact suggests that the behaviour of the theory near the deconfinement 
point is dominated by the macroscopic properties of the flux tube, which are largely 
independent of the precise ultraviolet details of the models and only relies on the in- 
frared effective string action. As a consequence some "super-universal" behaviours, 
like those that we have observed, are expected for dimensionless ratios of physical 
observables. Assuming this point of view, we expect that the effective models and 
ideas that we shall discuss in the next section should be valid for any gauge theory 
in 2+1 dimension. 

Let us conclude this section by noticing that in (3+1) dimensions, there is another 
interesting situation (which we will not discuss in this paper) where finite size effects of 
the type described in this paper can be observed. It is the cylindric geometry described 
in ref.s [||, || in which two spatial directions are chosen to be approximately of the same 
size as the inverse deconfinement temperature. 

This paper is organized as follows: in sect. 2 we give a short discussion of the flux tube 
model, in sect. 3 we describe the details of the Monte Carlo simulation and in sect. 4 we 
compare theoretical predictions and numerical results. 



2 The fluctuating flux tube model 

The flux tube model is based on the idea that in the confined phase of a gauge theory the 
quark-antiquark pair should be joined together by a thin, fluctuating flux tube 0. The 
simplest version of the model (which should be a good description when heavy quark- 
antiquark pairs are studied at large interquark separation R) assumes that the chromo- 
electric flux is confined inside a tube of small but nonzero thickness L c , that L c is constant 
along the tube (neglecting boundary effects near the quarks) and independent of the 
interquark distance. 

An immediate consequence of this picture is linear confinement: the interquark po- 
tential V(R) rises linearly according to the law V{R) = crR. 

A second, important consequence is that the string tension a and the effective cross- 
section of the flux tube A t are related by the law 

a= A t ' (1) 

where A t is defined as 

A _ UdAE) 2 



^~ J dA E 2 ' (2) 

In some cases the constant ct and, in particular, its dependence on the quark represen- 
tation, can also be determined (see ref. ||] for a discussion on this point and for further 
references). If the chromoelectric flux E is constant across the tube, then the simple 
dimensional relation At oc 2 ^ holds, where (d — 2) is the number of transverse 

dimensions. 
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Eq. (|l]) can be tested, for instance, by measuring the string tension on asymmetric 
lattices, with one space direction, say L s , smaller than the others. If L s > L c no effect on 
the string tension is expected. On the contrary, when L s < L c the flux tube is squeezed 
and the string tension increases. In particular, if we assume that the flux density is 
uniform inside the flux tube, then eq. (Q) suggests that the rising of the string tension is 
linear and obeys the law 

a(L s ) = ^a(oo) (3) 

where <r(oo) denotes the string tension in the uncompressed situation, namely for L s >> 
L c (in the following we will denote <t(oo) with a for brevity). 

This is the simplest possible assumption on the behaviour of the chromoelectric flux E 
across the tube and as a consequence the linear behaviour of eq.([|) is not at all mandatory 
and will have in general to be corrected when regions deep inside the flux tube are probed. 
Nevertheless, it has has been neatly observed recently in the case of the (2+1) dimensional 
SU(2) model The interesting aspect of such a linear behaviour implied by this uniform 
flux distribution is that it allows a rather precise determination of the flux tube thickness 
L c , which turns out to be in good agreement with other independent evaluations. For 
instance, in the above mentioned case of the (2+1) dimensional SU(2) gauge model, the 
value quoted by Teper j| agrees with the corresponding value obtained by Trottier and 
Woloshyn ||] with a different method. 

The fact that the flux tube has a finite non-zero thickness must be carefully taken into 
account to avoid systematic errors in the evaluation of the string tension from numerical 
simulations. It is in fact clear from the above discussion that only Wilson loops of size 
greater than L c must be used to extract the string tension in order to avoid the appearance 
of unphysical effects, in particular an artificial enhancement of the string tension. This is 
particularly important since small size loops, due to their small relative errors, dominate 
the fits and can strongly influence the determination of a. An important cross-check of 
these systematic errors is given by looking at the variation of x 2 /d.o.f. (from now on 
denoted by x 2 ) as a function of the minimal size L m i n of the Wilson loops considered 
in the fits. Only when L m i n is of the order of L c good Xr are obtained and stable values 
of the string tension are found . Notice also that these systematic deviations from the 
area law at small distances cannot be considered as lattice artifacts, because of their good 
scaling behaviour ||]. They could in principle be used as an independent way to estimate 
L c . 

Let us finally remark that the enhancement of the string tension at short distances is 
a general phenomenon: it can be precisely seen also in the case of 4d SU(2) and SU(3) 
gauge theories (see for instance fig. 2 of ref. [pj) and also when Polyakov line correlations 
are studied (see for instance fig.l and fig. 4 of ref. JTl|). 

A third, important, consequence emerges if one tries to take into account the quantum 
fluctuations of the flux tube. Below the roughening transition the Wilson loop follows the 
area law, which is the lattice counterpart of the "classical" linearly rising potential, but the 
continuum limit, hence the connection with the physical flux tube, can be achieved only 
after the roughening transition, where the quantum fluctuations of the surface bordered 
by the loop play a crucial role. It is by now generally accepted that these quantum 
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fluctuations can be effectively described by a massless two-dimensional free field theory. 
This idea traces back to the seminal work by Liischer, Symanzik and Weisz jl2| and has 
been discussed from then in several papers both in the context of lattice gauge theories 
(see for instance ref. fl3|] and references therein) and in the dual context of interfaces in 
3d spin systems (see for instance ref. []l4|] and references therein). Let us briefly review 
this approach and fix notations and conventions. 



2.1 Quantum fluctuations of the flux tube 

The starting point is the assumption that the fluctuations of the surface bordered by the 
Wilson loop are described by an effective Hamiltonian proportional to the change they 
produce in the area of the surface itself 



H, 



eff 



Li 



d,X2 



'1 + 



dh 



+ 



dh 
dx 2 



1 



(4) 



where the field h{x\,x 2 ) describes the surface displacement from the equilibrium position 
as a function of the longitudinal coordinates x\ and x 2 , L\ (L 2 ) is the size of the Wilson 
loop in the x\ (x 2 ) direction and a is the string tension []. 

The contribution to the Wilson loop expectation value due to surface fluctuations is 
then given by 

Z eff = tre- H «ff . (5) 

The Hamiltonian of eq. (Q) is too difficult to be handled exactly. However a crucial 
observation is that this theory can be expanded in the adimensional parameter (aL\L 2 ) 1 
and the leading order term is the gaussian model, which will be a good approximation 
when large enough Wilson loops are studied. Then we replace eq. (Q) with the aL\L 2 — > 00 
limit H -> H G 



I JO 



rL2 

dx\ 1 dx2 
Jo 



dh 

dx\ 



+ 



dh 
dx 2 



(6) 



Within this approximation, the integration over h implied in eq. (|5|) can be done exactly. 
The integral is divergent but can be regularized using, for instance, a suitable generaliza- 
tion of the Riemann ("-function regularization (see e.g. ]l5|]). The result depends only on 
the geometrical properties of the boundary. In particular, for a Wilson loop with fixed 
( "Dirichlet" ) boundary conditions along the loop, the gaussian contribution turns out to 
be 

Z G (L 1 ,L 2 )= ; C _ , T = i^ , (7) 



L 2 



where r/ denotes the Dedekind eta function 



r){r) 



I 24 



n=l 



(8) 



1 Note, as a side remark, that eq. (^) coincides with the Nambu-Goto string action in a special frame 
where only a subset of surface configurations are allowed. However the results of this section are also valid 
in a more general frame where all the possible configurations are taken into account. 
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c is an undetermined constant and we have assumed L\ > L 2 , without loss of generality. 

Thus, taking into account the area term aL\L 2 (coming from the classical or zeroth 
order contribution of the saddle point) and the undetermined perimeter contribution 
p{L\ + L 2 ) (which depends on the lattice regularization), the Wilson loop expectation 
value W(Li, L 2 ) may be written in the form 

W(L u L 2 ) = e- aLlL2+p(Ll+L2) Z G , (9) 

where the contribution Zq of the quantum fluctuations can be expanded as follows 



Z G {L X ,L 2 ) = q-*sJl +q + 2q 2 + 3q 3 + 5q± + ... (10) 
\/ L 2 

with q = exp(—2TTj^). Eq. ( flOj ) is symmetric under the exchange L\ <-> L 2 , as one can 
check by applying the modular transformation r — > — r _1 at eq. (^) (see e.g. Jl^]). 

As a consequence, a simple area-perimeter-constant law cannot fit the Wilson loop 
expectation values, and this shows up in high Xr' s an d unacceptable confidence levels^. 
On the contrary, if the contribution of the fluctuations is taken into account, impressive 
reductions of the Xr are found, and eventually acceptable confidence levels are reached 
(see for instance Q). This is up to now one of the strongest evidences for the correctness 
of this description of the fluctuations of the flux tube. 

One of the consequences of the above description in terms of a free bosonic field is 
that the mean width of the flux tube is expected to grow logarithmically as a function of 



the interquark distance [12], while there are reasons to believe that such a width should 
be constant (see for instance Ref . ]l6fl ) . Actually there is a slight modification of this 
theory which accounts for this fixed thickness, based on the observation that the free 
boson can be seen as a limit of a one-parameter family of conformal field theories, where 
this parameter can be simply related to the width Lq of the flux tube. 

In Ref. [16] it has been argued that the value of this parameter can be determined by 
matching the boundary conditions of the conformal theory with those of the underlying 
gauge field theory, combined with the obvious assumption that the flux tube cannot self- 
overlap freely. The latter assumption can be formulated more precisely in the Ising gauge 
theory by saying that the flux tube must sweep in its time evolution self-avoiding surfaces 



(see Ref. [13] for a discussion of this point). 



In this way one finds [16] Lq = yJn/Aa for a flux tube at zero physical temperature. 
In the following we shall not use this result, because we want to evaluate such a width 
by applying a completely independent argument to a different physical situation (the flux 
tube at the critical temperature T c ). We anticipate that the resulting value of Lq at T c 
turns out to be very near the one quoted above. 

Another consequence of this mentioned modification of the free boson theory is that 
the contribution due to the quantum fluctuations of the flux tube has now a different 
form: 



2 Notice, for completeness, that in the case of gauge theories with continuous gauge groups, perturbative 
contributions (for instance one gluon exchange terms) are also expected for small interquark separations. 
However these are completely independent from the above discussed infrared contributions and should 
be separately taken into account in fitting Wilson loop expectation values jlo|. Notice also that the 
above discussion requires non-trivial modifications if "fuzzy" or "blocked" versions of the Wilson loop are 
studied. 
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Z SA (Lx,L 2 ) = cq + ^ +ql + q 2 + . . . (11) 

where c is an undetermined constant and q = exp(—2irj^). 

Notice however that, while this modification has important consequences in the low 
temperature, confining regime, it is actually irrelevant in the high temperature phase we 
are discussing in this paper. Hence we will only refer to this choice when comparing our 
results with data taken in the low temperature phase. 



2.2 Finite Temperature 

The fact that the flux tube has a non-zero thickness becomes particularly relevant when 
the Wilson loop is studied at a finite temperature T, because this introduces into the 
game a new scale i.e. the length of the lattice in the time direction Lt = 1/T. For high 
enough temperatures, this length eventually becomes comparable with the finite thickness 
of the flux tube and the free field picture described above breaks down. 

Notice however that, since we are studying Wilson loops orthogonal to the time di- 
rection, this "temperature" interpretation is not mandatory. A space-like Wilson loop 
at high temperature is completely equivalent to an ordinary Wilson loop in a zero tem- 
perature environment, with the lattice size in one of the remaining space-like directions 
smaller than the other. This is exactly the situation described at the beginning of this 
section and studied by Teper in all the results listed there still hold with the simple 
exchange of L s with Lt. Exactly as before, choosing higher temperatures, namely looking 
at smaller and smaller values of Lj, we are actually probing the interior of the flux tube. 

In finite temperature LGT's the interquark potential can be extracted only by looking 
at the correlations of Polyakov loops. The surface bordered by two Polyakov loops under- 
goes a roughening transition exactly as in the Wilson loop case and the same arguments 
and techniques described above for the quantum fluctuations of Wilson loops apply also 



in this case [17, 18]. Using these results Olesen was able to predict the ratio between the 



deconfinement temperature and the square root of the zero-temperature (namely Lt = oo) 
string tension cr [18] 

T c / 3 

7^ = V^(d^2) (12) 

where (d — 2) is the number of transverse dimensions of the flux tube. This prediction 
turns out to be in good agreement with the results obtained with Monte Carlo simulations 
in (3+1) dimensions, both for the SU(2) and SU(3) gauge theories. The comparison with 
the numerical data for models in (2+1) dimensions is less good but the disagreement is 
anyway contained within 20%. For a review of these data see tab.III and ref. @. 



2.3 Vanishing of the quantum fluctuations 

Let us make at this point our main assumption. We assume that the quantum corrections 
Zq should disappear when the flux tube fills the whole lattice, namely when T > T c or 
L s < L c , depending on the geometry we are interested in. Intuitively this is equivalent to 
assume some sort of self-avoiding behaviour of the flux tube, since in this case , when the 
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flux tube fills the whole lattice there is no more space left for it to fluctuate. We will show 
below that this assumption is satisfied in the case of the (2+1) Z 2 model. This vanishing 
of quantum fluctuations can be described in a more rigorous way by noticing that the 
compactification in one lattice direction (say, L s ) naturally induces a compactification 
of the field h(xi,x 2 ) on a circle of radius R = The quantum field theory of a two- 
dimensional bosonic field compactified on a circle is by now rather well understood. In 
particular the spectrum of states and the partition function are exactly known (at least 
for the so called "rational models" for which R 2 is a rational number). It is thus possible 
to follow the behaviour of the quantum corrections as a function of R: Zq = Zq[K). 
What is interesting is that there are precisely four values of R (but only two of them 
are independent, the others being related by duality) for which the contribution of such 
quantum corrections vanishes (they correspond to the so called "topological field theories" , 
see ref. [13] for notations and bibliography on this subject). According to ref. [13, 
we can then predict that the value (let us call it Lq = 2ttRq) of the lattice size which 
corresponds to a zero-contribution point must be related to the (zero temperature) string 
tension as follows 

(13) 

Following our assumption we can thus say that moving toward higher temperatures 
corresponds in the flux tube effective model to a flow toward one of these zero-contribution 
points and that Lq must coincide with L c , the flux tube thickness. By comparing eq.s 
( |i~3| ) and (^) we then see that in (2+1) dimension the flux tube thickness L c coincides 
with the inverse deconfinement temperature. 

Eq.(13), and the fact that L c = 1/T C are in rather good agreement (within at most 
a 20% of deviation depending on the model) with the result of the simulations of the 
(2+1) dimensional SU(2) model ||, [|] and with our simulations of the (2+1) Z 2 gauge 
model (see tab. III). According to the subsection [d] of the introduction we suggest that 
this similarity should hold for any gauge theory in (2+1) dimensions. 

This similarity allows us to identify the deconfined regime with the compressed flux 
tube regime and to interpret the increasing of the space-like string tension in the deconfined 
phase as a signature of the compression of the flux tube. 



3 (2+1) Z2 gauge Monte Carlo simulation 

We compared our predictions with a set of high statistics simulations of the Z 2 gauge 
model in (2+1) dimensions. Let us briefly list some of the reasons for this choice (a more 
detailed discussion can be found for instance in ref.s H, [TT||), 

a] Due to the fact that this model is the dual of the 3d Ising model, very precise values 
of all bulk quantities (critical couplings and indices) are known. Moreover, in the 
region that we study an excellent agreement with the scaling laws has been found. 

Choosing the standard normalization for the action 

S = -p'£U p > (14) 
v 
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where U p denotes the product of Z2 elements associated to the links belonging to the 
plaquette p, we have a roughening transition at /3r ~ 0.4964, a (zero temperature) 
deconfinement transition at f3 c ~ 0.7614 and a critical index for the correlation 
function v ~ 0.63. Precise information on the finite temperature deconfinement 
transition can be found, for instance, in pi] ]. The scaling region of this transition 
starts approximatively at Lt = 4. Data obtained with the Monte Carlo renormal- 
ization group approach for L t = 4 and L t = 8 give: (3 c {L t = 4) = 0.7315(5) and 
j3 c {L t = 8) = 0.7516(5) |H] (where the notation f3 c {L t ) denotes the value of at 
which deconfinement occurs in a lattice with temporal size Lt). Defining the critical 
temperature T c as 

' T§ (f3 c -f](L t )y , (15) 



we have Tfi = 2.3 (1). 



b] The model is simple enough to allow high statistics simulations, but shares the same 

infrared behaviour (in the (3 < (5 C region) with non abelian lattice gauge theories. 
In particular, near the finite temperature deconfinement point, the Svetitsky-Yaffe 
conjecture tells us that Z2 and SU(2) models should be in the same universality 
class. 

c] Recently the interface behaviour in 3d statistical models has attracted a lot of interest. 

In d=3 the physics of interfaces is exactly equivalent, through duality, to the physics 
of Wilson loops and, in the particular case of the Ising model, not only theoretical 
tools but also numerical results can be borrowed from one context to the other (see 
ref. [|14|] for a discussion of this problem). In particular, high precision estimates of 
the string tension a (equivalent to the interface tension) exist. In the scaling region 
the string tension behaves as 

a(J3) = a (/3 C - Hf v • (16) 

It turns out that this law is very well fitted by the existing data [14] with cjq = 
3.70 (4). By inserting the value for the critical temperature obtained from eq. (|l5|), 
we can construct the adimensional ratio: Tq/^/oq = 1.19 (6). 



3.1 The simulation 

The simulation was performed on a 48 2 x Lt lattice {Lt = lattice spacing in the time 
direction) with periodic boundary conditions in all space-time directions. A standard 
heat-bath algorithm was used to update links. Four values of /3, with various choices of 
Lj, were studied. They are listed in tab. I . 



p 


L t 


Lc(P) 


0.7420 


13 


5.2 (2) 


0.7500 


10,12 


7.3 (3) 


0.7525 


2-8 


8.5 (4) 


0.7585 


2-16 


17.2 (7) 
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Tab. I. Set of measured data. In the third column, the corresponding inverse critical 
temperatures (in unities of the lattice spacing) are presented. 

The reason for these choices is the following. The data at (3 = 0.7585 have been 
measured to make quantitative estimates of the temperature dependence of a(Lt). Since 
in this case the finite temperature deconfinement point is at Lt ~ 17, a wide range of 
values of Lt can be used to check the predicted behaviour of the string tension. The data 
at (3 = 0.7525 have been taken to check the scaling behaviour and the overall reliability 
of our results. In both cases all deconfined values of L t were studied, namely 2 < L t < 8 
for = 0.7525 and 2 < L t < 16 for = 0.7585. 

Finally, to make a comparison with the above data, we studied three samples in 
the confining phase at intermediate finite temperatures (from T/T c ~ 0.4 to T/T c ~ 
0.73): these have been chosen in order to have low enough temperature to guarantee an 
uncompressed behaviour of the flux tube (in particular it can fluctuate freely). The values 
of (3 ([3 = 0.742 and (3 = 0.750) have been tuned so as to have values of the string tension 
comparable with those of the compressed case. In particular, a((3 = 0.742) ~ a(f3 = 
0.7525, L t = 4) ~ a{(3 = 0.7575, L t = 4) and a{(3 = 0.7575, L t = 6) < a{(3 = 0.750) ~ 
a((3 = 0.7525, L t = 6) < a((3 = 0.7575, L t = 5), (see tab.IIa,6,c). 

All Wilson loops orthogonal to the time direction W(L\, L2) in the range 2 < L%, L2 < 
20 were measured. Also this choice requires some explanation. It is well known that 
improved versions of the Wilson loop operator [ p^ ] give more precise results, but notice 
that only in the case of ordinary Wilson loops the gaussian determinants described in 
the previous section can be evaluated exactly. Since our goal is to show the fate of the 
quantum fluctuations of the flux tube across the deconfinement transition, more than to 
have precise evaluations of, say, L c , we decided to concentrate on ordinary, not improved, 
Wilson loop expectation values. 

3.2 The cross-correlation problem 

The major problems one has to face for these values of f3 are the critical slowing down 
and the huge cross-correlations among Wilson loops. The first problem was kept under 
control separating each measurement on the lattice with 32 sweeps. The second was taken 
into account, as usual, by weighting the data in the fitting procedure with the inverted 
cross-correlation matrix . In the most severe cases, namely those near the critical point 
(L t = 6 — 8 for (3 = 0.7525, and L t = 13 — 16 for (3 = 0.7585) a scattering procedure was 
also implemented, measuring in each iteration only one single Wilson loop of fixed size 
and scattering the measure of the others in the Monte Carlo time0. Let us stress that all 
these steps were absolutely crucial to obtain reliable confidence levels in the fits. Moreover 
in these last cases further runs with lattice size 96 2 x Lt were performed, measuring loops 
up to 40 x 40, to check the reliability of our results. Notice finally, as a side remark, that 
improved Wilson loop estimators experience even worse cross-correlations. 
All the quoted errors were obtained with a standard jacknife procedure. 

3 In these cases the cross-correlation matrix was so flat, within the statistical errors, that the reliability of 
the whole inversion procedure and cross-correlated fit calculation, without a proper scattering procedure, 
was rather doubtful. 
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4 Results and conclusions 



For each value of j3 and L t we fitted the Wilson loop expectation values with a pure area 
law 

W(L 1 ,L 2 ) = exp{-aL 1 L 2 + p{L l + L 2 ) + c} (17) 

(we shall call this choice "type 1" fit in the following) and with an area law corrected in 
order to take into account the quantum fluctuations of the Wilson loop surface 

W(L 1 ,L 2 ) = Z G (L 1 ,L 2 ) exp{-uL 1 L 2 + V [L X + L 2 ) + c} (18) 

("type 2" fit in the following). For the three samples in the confined phase we fitted the 
data also with the modified version of the gaussian contribution proposed in eq. ( |Tl|) 

W(Lx,L 2 ) = Z SA (Lx,L 2 ) ex V {-oL x L 2 + V {L X + L 2 ) + c} (19) 

("type 3" fit in the following). 

We performed the fits setting a lower threshold Rt in the size of the Wilson loops. To 
be precise we constructed the following subsets of our Wilson loops samples 

S{R t ) = {W{L U L 2 ), L X >L 2 > R t } (20) 

choosing Rt = Lt in the deconfined phase. In the confined region we chose instead 
Rt = L c ((3) (namely the flux tube thickness, extrapolated through scaling at those values 
of (3). This means R t = 4 for (3 = 0.742 and R t = 6 for = 0.750. Notice that, as a 
consequence of this cutoff procedure, severe problems of precision exist near the critical 
temperature. The results of the fits are collected in tab.IIa, 6, c and fig. 1,2. 
Let us make some comments on these data. 

Xr behaviour 

In the deconfined phase type 1 fits have in general better confidence levels than type 
2 fits (see tab.IIa and tab. 116). This is quite evident in the high temperature regime 
T > 2T C , while near the critical point the two Xr' s almost coincide, indicating that 
the quantum fluctuations corrections in these region are smaller than the precision 
of our data and could not be detected in any case even if they are present. This 
behaviour is even more distinct if compared with analogous fits for the two samples 
of data at low temperature, in the confined region (see tab. lie), where definitely 
better confidence levels are obtained with fits of type 2 and 3. This is a first 
signature of the vanishing of quantum fluctuations. 

T c <T < 2T C 

Let us concentrate on the data at (3 = 0.7585 (tab.IIa). We fitted the values of 
u(Lt) in the range 9 < Lt < 16 according to the linear law of eq. (||). As it can 
be seen in fig.l, the linear law is in good agreement with the data, and in fact the 
fit shows a good confidence level: C.L.=98%. As a result we obtain L c = 18.5 (9), 
which implies: L c = 1.07(5)(5)T C _1 and L c ^/oq = 0.90(5)(2), where the second error 
in these relations is due respectively to the uncertainty in the critical temperature 
and the string tension. 
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T/T c 


a 


O"0 


xi 


xi 


2 


8.6 (4) 


0.120 (3) 


189. (5) 


1.3 


7.1 


3 


5.7 (3) 


0.0480 (2) 


75.6 (3) 


0.71 


8.1 


4 


4.3 (2) 


0.0244 (3) 


38.8 (5) 


0.85 


5.1 


5 


3.45 (17) 


0.0150 (4) 


23.6 (6) 


1.2 


1.8 


6 


2.87 (14) 


0.0100 (3) 


15.8 (5) 


0.91 


1.8 


7 


2.46 (12) 


0.0083 (5) 


13.1 (8) 


1.3 


1.3 


8 


2.15 (11) 


0.0068 (4) 


10.7 (6) 


0.78 


1.3 


9 


1.92 (10) 


0.0049 (4) 


7.7 (6) 


0.90 


1.08 


10 


1.73 (9) 


0.0044 (4) 


6.9 (6) 


1.1 


1.4 


11 


1.57 (8) 


0.0037 (4) 


5.8 (6) 


0.72 


0.80 


12 


1.44 (7) 


0.0038 (4) 


6.0 (6) 


1.01 


1.06 


13 


1.33 (7) 


0.0033 (5) 


5.2 (8) 


0.90 


0.90 


14 


1.23 (6) 


0.0033 (5) 


5.2 (8) 


0.52 


0.74 


15 


1.15 (6) 


0.0025 (8) 


3.9 (1.3) 


0.46 


0.48 


16 


1.08 (5) 


0.0022 (9) 


3.5 (1.4) 


0.85 


0.85 



a 





T/T c 


a 


0"o 


xi 


xi 


2 


4.2 (2) 


0.128 (2) 


49.1 (8) 


0.79 


2.73 


3 


2.84 (14) 


0.0520 (4) 


19.9 (2) 


0.81 


1.68 


4 


2.13 (10) 


0.0270 (10) 


10.0 (4) 


1.08 


1.44 


5 


1.70 (8) 


0.0172 (5) 


6.6 (2) 


0.96 


1.00 


6 


1.42 (7) 


0.0137 (7) 


5.3 (3) 


0.70 


0.81 


7 


1.22 (6) 


0.0112 (7) 


4.3 (3) 


0.80 


0.80 


8 


1.06 (5) 


0.0085 (9) 


3.3 (3) 


0.70 


0.70 



b 



p 


L t 


T/T c 


a 


O"0 


xi 


xi 


xi 


0.742 


13 


0.40 (2) 


0.027 (2) 


3.88 (20) 


1.7 


0.85 


0.80 


0.750 


12 


0.61 (3) 


0.0132 (10) 


3.71 (29) 


1.4 


1.1 


1.02 


0.750 


10 


0.73 (3) 


0.0127 (11) 


3.56 (32) 


1.4 


1.1 


0.89 



c 



Tab. II. (a,b,c) Values of the string tension extracted from Wilson loops at 
P = 0.7585 (a); (3 = 0.7525 (b); (3 = 0.7420, 0.7500 (c). In the first two columns 
the lattice size in the time direction and the corresponding temperature (in units of 
the critical temperature) are reported. The last columns contain the reduced x 2 of 

the fits of type 1 and 2 (in tab. lie also fits of type 3 are considered). The third 
column contains the string tension cr(Lt), extracted from fits of type 1 in the cases 
a,b and fits of type 3 in the case c. In the fourth column are reported the 
corresponding scaling values o~o(Lt), according to eq. 
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A slightly lower confidence level (35%) is obtained if the same fit is performed 
on the (3 = 0.7525 data (tab. 116), with L c = 8.7(5), which corresponds to L c = 
1.02(5)(5)T- 1 and L c ^ = 0.86(5)(2). 

These results are summarized in tab. Ill, where they are also compared with the 
corresponding values for the SU(2) model Q || and with our predictions eq.s(|l2]) 
and (0). 



model 






L C T C 


SU(2) 


~ 1.09 


1.12(1) 


~ 1.22 


z 2 


0.90(7) 


1.19(6) 


1.07(10) 


FTM 


1.024 


0.97 


1 



Tab. III. Comparison between Z 2 and SU{2) gauge theories in (2+1) dimensions. 
In the last row are reported the flux tube model (FTM) predictions eq.s^!2\) and 

m 



T > 2T C 

In this region the linear behaviour of eq. (||) is completely lost. This indicates the 
appearance of some non-trivial structure deep inside the flux tube. The data show 
a good scaling behaviour as a function of T and suggest that the only remaining 
physical scale in this regime is the temperature. 

This can be seen by looking at fig. 2. Since the errors in this region are very small, 
the Xr becomes a very efficient tool to select among various possible behaviours. 
Indeed, as can be seen in fig. 2, both the data at (3 = 0.7525 and those at (3 = 0.7585 
have a very precise power law behaviour 



a(L t ) = aL? 



(21) 



and the two values of a and a at the two /3's are almost compatible within the 
errors. 

The output of the fits is 

(3 = 0.7525, range: 4 > L t > 2, C.L.= 66%, a = 0.58(2), a = -2.19 (4) , 
(3 = 0.7585, range: 6 > L t > 2, C.L.=75%, a = 0.59(2), a = -2.28 (4) . 

4.1 Vanishing of quantum fluctuations 

In order to have an independent check of the vanishing of quantum fluctuations in the 
deconfined phase we constructed the ratios 



C{L,n) 



W(L,L) 



(22) 



W(L + n,L-n) ' 

If the Wilson loops are described by a pure area law then C(L, n) does not depend on L 

C(L,n)=e- an2 (23) 
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while, if a contribution coming from the flux tube fluctuations is present, a decreasing 
function of L is expected 



Note that, being only two Wilson loops involved in these ratios, one has much smaller 
errors than, for instance, in the case of Creutz ratios. We decided to compare samples 
of data with (almost) the same string tension and the same expected flux tube thickness 
in order to minimize systematic errors (due, for instance, to different lower thresholds 
in fitting the data). Hence we compared the Wilson loop expectation values at j3 = 
0.7585, L t = 4 and = 0.7525, L t = 4 with those at j3 = 0.7420, L t = 13. It can be 
neatly seen by looking at fig. 3 that the two samples in the deconfined phase (compressed 
flux tube) have almost coinciding ratios and show no contribution coming from flux tube 
fluctuations, while the ratios of the sample in the confining region lie in between the 
slopes predicted by eq. (|3|) an d eq. (p5|). The same picture is confirmed by fig.4, where 
the samples at f3 = 0.7585, L t = 5,6 and (3 = 0.7525, L t = 6 are compared with those 
at (5 = 0.7500, Lt = 10, 12. Notice that in this case the lower threshold is Rt = 6, 
nevertheless also the data at lower values of R have been included in the figure, in order 
to show the enhancement of the string tension below the threshold, as discussed in sect. 2. 
Note also that the proposal of eq. ( |IT| ) and eq. ( |25| ) for the flux tube fluctuations seems to 
describe better the data exactly in that region in which the self-avoiding constraint plays 
a major role. 

Thus we can conclude that the information extracted from C(L,n) is in complete 
agreement with those coming from the Xr °f the fits and support the picture of the 
vanishing of quantum flux tube fluctuations in the high temperature phase. 

4.2 Comparison with the SU(2) gauge theory 

A remarkable feature of our results is their similarity with those obtained by Teper in the 
case of SU(2) in d=2+l dimensions (see tab. III). Indeed, this also is a consequence of the 
flux tube picture, which, at least in a first order approximation, does not depend on the 
ultraviolet details of the model (not even on the fact that the gauge group is discrete or 
continuous, abelian or non-abelian), but only on its infrared behaviour (hence essentially 
on the existence of a confining phase beyond the roughening transition and on the space 
time dimensions). Let us stress that the comparison made in tab. Ill is only a partial 
check of this supposed universality, since both the SU(2) and the Ising model have the 
same center. It would be quite interesting to check it also in other cases, and in particular 
for the SU(3) model. 
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(24) 



or, following eq. (|IT|) 




(25) 
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Figure Captions 



Fig.l. The string tension in scaled units (o~o(Lt) in tab.IIa ,b,c) is plotted as 
a function ofT/T c . Data corresponding to tab.IIa (triangles), tab. lib (squares), tab. lie 
(crosses), are reported. The horizontal line is the uncompressed value of the string tension 
Co = 3.70. The other line is the best fit to eq. (0) for the sample of data at (3 = 0.7585 in 
the range 9 < L t < 16. 

Fig. 2. The string tension (cr(Lt) in tab.IIa ,b) is plotted as a function ofT (1/Lt in 
tab.IIa,b). Data corresponding to tab.IIa (triangles) and tab. lib (squares), are reported. 
Both cr(Lt) and T are in log. scale. The continuum line is the best fit to eq. for the 
sample of data at (3 = 0.7585 in the range 2 < Lt < 6 . 

Fig. 3. The logarithm of the Wilson loop ratios of eq. (fli|), with n=l, are plotted 
for 13 = 0.7420, L t = 13 (triangles), (3 = 0.7525,L t = 4 (squares) and (3 = 0.7585, L t = 
4 (crosses). The corresponding value of the string tension (as it is listed in tab. II) is 
subtracted for each sample. The continuous line is the expected result if a pure gaussian 
correction is taken eq. §fi\), the dashed line is the self-avoiding case eq. (p?3[). 

Fig. 4. Same as fig. 3. The ratios are evaluated at step n=2. The various samples plot- 
ted are: (3 = 0.7500,L t = 10 (diamonds), (3 = 0.7500,L t = 12 (circles), (3 = 0.7525,L t = 6 
(squares), (3 = 0.7585, Lt = 5 (crosses), (3 = 0.7585, Lt = 6 (triangles). 
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